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1 . Introduction 

This chaptei0 deals with the following basic questions, which sometimes surprise or even 
irritate the reader, namely: Which are the intrinsic properties of the graphene layers inside 
the graphite structure? Are their transport properties quasi-two dimensional and are they 
comparable with those of single layer graphene? Which is the coupling between ideal, i.e. 
defect free, graphene layers inside graphite? Which is the influence of defects and interfaces 
within the graphene and graphite structure? These basic questions remain in the literature 
still unanswered to a large extent mainly because the influence on the transport properties of 
defects and sample or crystal sizes was not systematically taken into account. 
A large amount of the published interpretations on the experimental transport data of real 
graphite samples relied on the assumption that these were intrinsic properties. Let us 
start with a rather simple example. In the last fifty years scientists flooded the literature 
with reports on different kinds of electronic measurements on graphite samples, providing 
evidence for carrier (electron plus hole) densities per graphene layer in grap hit and at l ow 
temperatures hq > 10^" cmr^_as one c an read in the standard book from Kellv (1981) or 
the old publication from iMcClurd il964h . or from more recent work by Grun eis et al. 1 20081) 
(wo - 10^^ cm"2) or lKimiar et alj bOlOh (wn ~ 2.4 x lO" cm"^) ^j^g concentration 

is not a constant and does apparently depend on the measured sample? The reader may ask 
then why at all is the knowledge of the intrinsic carrier concentration ng so important? 
The carrier concentration is one of those basic parameters needed to estimate several others 
necessary to get a reliable picture of the electronic band structure and from this to understand 
all intrinsic (or not) transport properties. For example, two-dimensional (2D) calculations 



^ This work is supported by the DFG under DFG ES 86/16-1. A.B. is supported by the ESF 
"Nachwuchsforschergruppe" "Funktionale multiskalige Strukturen" and S.D. by the Graduate school 
of Natural Sciences of the University of Leipzig "BuildMona". 

^ In this chapter and to facilitate a direct comparison between the carrier density of a single graphene 
layer and that one inside the graphite structure, the carrier density per area will be used. This can be 
easily obtained multiplying n, the three dimensional carrier density, by the distance between graphene 
layers in graphite flo ~ 0.335 nm. This estimate assume that the coupling between graphene layers 
inside graphite is very weak and that the electron transport is mainly on the two dimensional graphene 
layers. 
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assuming a coupling 70 between nearest in-plane neighbors C-atoms on the graphene plane 
give a carrier density (per C-atom) w(T) = (0.3 . . . 0.4) (^8^/70)^ (7o — 3 eV and T is the 
temperature) jKellyl Il98lh . Because aU experimental values obtained from bulk graphite 
samples indicated a finite n{T ^ 0) = mq > then the straightforward and "easiest" solution 
to solve this "disagreement" is to start including more free parameters in the tight-binding 
electronic band structure calculations. 

For example, introducing a new coupling 71 between C-atoms of type a in adjacent planes one 
obtains n{T) = a{yi / Jq)T + b{T / jq)'^ + c{T^ / j^ji) + . . .{a,b,c,. . . are numerical constants), 
where the "accepted" value for 71 ~ 0.3 eV. Also in this case n(T ^ 0) — )• 0. We stress that 
neither in single layer graphene nor in graphite such T— dependence was ever reported 0, 
i.e. a large density backgroimd «o was always measured and assumed as "intrinsic" without 
taking care of any influence from lattice defects or impurities . To fit exp eriment al data and 
obtain a finite Fermi energy Ef - in the simplest case Ep ~ 72 jPillon eta l.. 1977; Kellyl[l98lh 
- up to seven free parameters were and still are introduced, even for carrier density as small 
as w ~ I — 8 X 10' I cm~ ^(EF ~ —29 meV) as o btained recently from magnetotransport data in 
bulk pyrolytic graphite jSchneider et al J. [20091) . 

Taking into acc ount the exhaustiv e experience accumulated in gapless or narrow gap 
semiconductors jTsidilkovskii Il997h we should actually expect that the measured carrier 
density wq ^ 10^^ cm^^ is not intrinsic of the graphite structure but it is influenced by 
impurities and/ or defects in the graphite /graphene lattice. The reader should keep in mind 
that a carrier density of the order of 10* cm^^ means one carrier in 1 /im^ graphene area, which 
could be produced actually by having a single vacancy or impurity in the same graphene area, 
in case one carrier is generated by one of these defects, as experimental ( Arndt et al., 2009) and 
theoretical iStauber et alJ.l200^ work suggests. Experimental evidence published recently and 
partially reviewed in this chapter speaks against an intrinsic origin of - even a part of - the 
measured mq in graphite samples, casting doubts on the relevance of related electronic band 
structure parameters obtained in the past. On the other hand this new knowledge will help 
significantly to clarify observed transport phenomena. 

In Section|3]of this chapter we describe a method that one can use to obtain the carrier density, 
the mean free path and the mobility of the carriers inside graphite without free parameters. In 
that Section we review a systematic study of the transport in small multigraphen^ samples 
that reveals room-temperature mobility values (~ 6 x lO'' cm^V^^s^^) per graphene layer 
inside graphite, which overwhelm those reported in literature for single graphene layers, 
indicating the higher quality of the graphene layers inside ideal graphite. This quality is also 
reflected by the extremely low room-temperature carrier density (~ 7 x 10* cm^^) obtained 
for good, but certainly not ideal, quality multigraphene samples. These studies indicate that 
ballistic transport with mean free path in the micrometer range is possible in graphite at room 
temperature. 

In Section |2] we describe the main characteristics of bulk and multigraphene samples, 
their characterization using transmission electron microscopy (TEM), electron backscattering 



^ It is interesting to note that the carrier concentration obtained in bulk graphite bv lGarcia et alj j2008l) , 
using an original and parameter-free method to determine it and the mean free path, can be fitted up to 
~ 200 K by n[cm"2] ^ uq + lO^T^ + 7.5 x lO^T^ with T in [K] and ng ~2x 10* cm"^ y^e same data, 
however, can be also well explained by a semiconducting-like exponential function exp(— E^/2T) with 
an energy gap Eg ~ 50 meV. 

* We use this word to refer to graphite samples of not more than a few micrometers in length and width 
and thickness below ~ 100 nm. The reason for this kind of geometrical restriction will become clear in 
Section|2] 
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diffraction (EBSD), electronic transport and Raman spectroscopy. We show the correlations 
between the internal microstructure and sample size - lateral as well as thickness from 
millimeter size graphite samples to mesoscopic ones, i.e. tens of nanometer thick 
multigraphene samples - and the temperature (T) and magnetic field (B) dependence of the 
longitudinal resistivity jo(T, B). Low energy transmission electron microscopy reveals that the 
original highly oriented pyrolytic graphite (HOPG) material - from which the multigraphene 
samples were obtained by exfoliation - is composed of a stack of ~ 50 nm thick and 
micro meter long crystalline regions separated by interfaces running parallel to the graphene 
planes iBarzola-Ouiquia et al .l 2008). We found a qualitative and quantitative change in the 
behavior of p(T, B) upon thickness of the samples, indicating that their internal microstructure 
is important. The overall results described in sections |2] and |4] indicate that the metallic-like 
behavior of p(T) at zero magnetic field measured for bulk graphite samples is not intrinsic of 
ideal graphite. 

The influence of internal interfaces on the transport properties of bulk graphite is 
described in detail in Section |4] of this chapter. We will show that in specially prepared 
multigraphene samples the transport properties show clear signs for the existence of granular 
superconducti vity within the graphite interfaces, which existence was firstly reported 
by iBarzola-O uiquia et al.l J2008h . Based on the results described in Section |4] we argue that 
the superconducting-insulator or metal-insulator transition (MIT) reported in literature for 
bulk graphite is not intrinsic of the graphite structure but it is due to the influence of these 
interfaces. 

2. Samples characteristics 
2.1 Sample preparation 

In order to systematically study the transport properties of ideal graphite and compare them 
with those of bulk graphite samples measured in the past, we need to perform measurements 
in different tens of nanometer thick multigraphene samples of several micrometer square area. 
The samples presented in this chapter were obtained from a highly oriented pyrolytic graphite 
(HOPG) bulk sample with a mosaicity (rocking curve width) of ~ 0.35° ±0.1° from Advanced 
Ceramics company. This material does not only guaranty high crystalline quality but also 
allows us to easily cleave it and obtain up to several himdreds of micrometers large flakes 
with thickness from a few to several tens of nanometers. The starting geometry of the bulk 
graphite material for the preparation of the flakes was ~ 1 mm^ and ~ 20 jim thickness. 
The selected piece was located between two substrates and carefully pressed and rubbed. As 
substrate we used p-doped Si with a 150 nm SiN layer on top. The usefulness of the SiN 
layer on the Si substrate is twofold: firstly the multigraphene flake on it shines with high 
contrast by illuminating it with white light allowing us to use optical microscopy to select the 
multigraphene samples. 

Immediately after the rubbing process we put the substrate containing the multigraphene 
films on it in a ultrasonic bath during 2 min using high concentrated acetone. This process 
cleans and helps to select only the good adhered multigraphene films on the substrate. After 
this process we used optical microscopy and later scanning electron microscopy (SEM) to 
select and mark the position of the films. Figure [T] shows four of the investigated samples of 
different micrometer length and tens of nanometer thickness. 

For the preparation of the electrical contacts we used conventional electron beam lithography. 
The contacts were done by thermal deposition of Pd (99.95%) or Pt/Au (99.5%/99.99%) 
in high vacuum conditions. We have used Pd or Pt/Au because these elements do not 
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Fig. 1. (a) to (d): Optical microscope pictures of four multigraphene samples prepared as 
described ir\ the text. The dimensions of the samples can be seen directly at the pictures. The 
thickness of the samples are: (a) 60 nm, (b) 55 nm, (c) 10 nm, and (d) 85 nm. 

show any significant Schottky barrier when used with graphite. This has been checked by 
I — V measurements in the temperature range of the measurements used in this work. The 
advantage of our preparation method lies in the easy way to do (one rubbing process is 
enough to produce samples, and from the initial flake is possible to produce ten of substrates 
containing multigraphene samples), it avoids contamination or surface doping avoiding the 
contact with materials as with the Scotch-tape method. For the transport measurements the 
substrate with the sample was glued on a chip carrier using GE 7031 varnish. The contacts 
from the chip carrier to the electrodes on the sample substrate were done using a 25 jim gold 
wire fixed with silver paste. 

2.2 Transmission electron microscopy and electron backscattering diffraction 

The scanning electron microscope (SEM) pictures, electron beam lithography, lamella 
preparation and electron backscattering diffraction (EBSD) of the investigated samples were 
done using a Nova NanoLab 200 dual beam microscope from the FEl company (Eindhoven). 
A HOPG lamella was prepared for transmission electron microscopy (TEM) using the in-situ 
lift out method of the microscope. The TEM lamella of HOPG was cut perpendicular to 
the graphene layers. The electron transmission parallel to the graphene layers provides 
information on the crystalline regions and their defective parts parallel to the graphene 
layers. We obtained thin lamellas of aroimd 200 nm thickness, ~ 15 //m length and ~ 5 }im 
width. After final thinning, the sample was fixed on a TEM grid. A solid-state scanning 
transmission electron microscopy (STEM) detector for high-resolution analysis (included in 
our microscope) was used. The voltage applied to the electron column was up to 30 kV and 
the currents used were between 38 pA to 140 pA. 

Figure|2]shows the bright (a) and dark field (b) details obtained with the low-voltage STEM at 
18 kV and (d) the bright field picture of a different lamella obtained at 30 kV. Figure|2l|c) shows 
a blow out of a detail of (a). The different gray colours indicate crystalline regions with slightly 
different orientations. The images indicate that the average thickness of the crystalline regions 
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Fig. 2. Transmission Electron Microscopy pictures taken parallel to the graphene layers of 
two HOPG lamella. The c-axis is perpendicular to the clearly observable stripes of different 
gray colours, each representing a crystalline region with a slightly different orientation. The 
arro ws in (a) and (b) indicate 400 nm length scale, in (c) 100 nm and in (d) 1 }tm. Adapted 
from lBarzola-Ouiquia et alj j2008h . 

is 60 ± 30 nm. One can clearly resolve the interfaces running perpendicular to the c-axis of the 
layers, parallel to the graphene layers. In Fig.|2tc) we also realize the end of crystalline regions 
along, normal or with a certain angle respect to the graphene layers direction. 
Other experimental technique used to characterize the crystal orientation, defects and 
grain boundaries is the electron backscattering diffraction (EBSD). The measurements were 
performed with a commercially available device (Ametek-TSL) included in our microscope. 
In this setup a cleaved bulk HOPG sample under investigation was illuminated by the SEM 
beam at an angle of 70° and the diffracted electrons were detected by a fluorescence screen 
and a digital camera. The included software was used to calculate the orientation of the 
crystalline regions within the HOPG surface as a function of the electron beam position. 
Figure |3] shows the grain distribution at the near surface region of a HOPG sample where 
the in-plane orientation is recognized by the (bluegreen) colour distribution. We recognize in 
this figure that the typical crystal size (on the (a,b) plane) in our HOPG samples is between a 
few fim and ~ 20 fim. Taking into account the TEM pictures shown in Fig.|2]and the EBSD 
one in Fig.|3]we conclude that single crystalline regions in a HOPG good quality sample is not 
more than ~ 20 iim long and less than ~ 60 nm thick. If the interfaces between the crystalline 
regions as well as the defects in the crystalline structure in graphite have some influence on the 
transport properties we then expect to see a change in the behavior of the transport properties 
between samples of thickness of the order or less than the average thickness of the crystalline 
regions and of planar size below ~ 10 fim. 
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Fig. 3. Electron backscattering diffraction pichires taken from the surface of a HOPG sample. 
The area scanned is :i 80 x 65 }im^ and the depth less than 10 nm. 

2.3 Influence of the internal microstructure on the temperature and magnetic field 
dependence of the longitudinal resistivity 

In order to investigate the intrinsic properties of graphite and taking into account the internal 
structure of bulk graphite samples, it appears obvious to study tens of nanometers thick 
graphite samples. iBarzola-Ouiquia et al.l feOOSh correlated the internal microstructure and 
sample size with the temperature and field dependence of the electrical resistivity. Whereas 
HOPG or graphite samples with thickness larger than ~ 50 nm show a metallic-like behavior 
in the electrical resistivity vs. temperature, in the case of tens nanometer thick graphite 
samples, it steadily increases decreasing temperature, see Fig.lHa). 




Temperature (K) Magnetic field (T) 



Fig. 4. (a) Normalized resistance versus temperature at zero applied fields for two 
multigraphene samples, namely: L5 (red dots) and L7 (empty black squares). The values of 
thickness and resistivity at 4 K are included in the inset table. The magnetoresistance 
(defined as shown in the y-axis) at three different temperatures are shown for samples L5 (b) 
and L7 (c). The magnetic fie ld was applied parallel to the c-axis. Adapted from 
iBarzola-Ouiquia et all i2003) . 
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The true temperature dependence of the resistivity in an ideal, defect-free multigraphene 
sample shou ld be semiconducting-like, as it is expected f or an ideal semimetal with zero 
or small gap. iDusari. Barzola-Ouiquia. Esquinazi & Garcial j201lh noted that the temperature 
dependence of several multigraphene samples shows a semiconducting-like behavior with a 
saturation or a weak maximum at low temperatures, this last behavior due to the contribution 
of interfaces, internal as well as with the substrate or at the free sample surface, parallel to 
the graphene layers, see also Section H) It becomes now clear that the contribution of the 
interfaces to the measured conductance of large samples can overwhelm the one coming 
from the intrinsic graphene layers in thick enough samples. Recently done studies on several 
graphite samples of different thickness indicate the existence of an int rinsic energy gap o f the 
order of 40 meV in the graphene layers within ideal Bernal graphite dGarda et al.l,l20lil) . We 
note that experiments on samples with carrier concentration < lO' cm^^ are hardly reported 
in literature. At such low-carrier densities, as it appears to manifest in graphite, electron 
correlations and possible localization effects should be considered. Electron interactions 
are large and for a small-enough carrier density, t he expected screenin g will be very weak 
promoting therefore the existence of an energy gap forut & Lahdel,l200^ . 
According to literature the structural quality of graphite samples can be partially resolved by 
investigating the Raman D-line at ~ 1350 cm^^. However, also the edges of the samples as 
well as the borderlines between regions of different thicknesses may contribute to the D-band 
signal. The Raman spectra between 1300 and 1700 cm~^ have been measu red at different 
positions of several multigraphene samples by iBarzola-Ouiquia et al.l j2008h . Those results 
show that in thin multigraphene samples with similar semiconducting-like behavior in the 
temperature dependence of the resistivity, different amplitudes of the Raman D-line peak are 
measured. This indicates that the disorder related to this Raman line does not appear to affect 
strongly the temperature dependence of the resistivity. 

As the absolute value of the electrical resistivity depends also on the measured sample 
thickness, see Section 12.41 it seems clear to assume that the metallic-like behavior is not 
intrinsic to ideal graphite but it is due to the influence of the interfaces inside the graphite 
samples with large enough thickness. Magnetoresistance (MR) results follow the behavior 
observed in the electrical resistivity. If a metallic behavior is obtained (thicker samples) then 
the MR decreases with increasing temperature and it shows SdH oscillations, see Fig.lHc). If 
the graphite sample shows a semiconducting-like behavior (thinner samples) then the MR is 
clearly smaller in magnitude as well as it decreases decreasing temperature for fields above 1 T, 
see Fig.lHb). This figure also shows the absence of the SdH oscillations when the thickness of 
the graphite sample is small enough. 

We note that Kohler's rule does not apply in multigraphene and HOPG samples and one of 
the reasons might be the contribution of the sample internal structure and interfaces. Defects 
and interfaces may influence the dimensionality of the transport and might be the origin of 
localized granular superconductivity. Results presented so far correspond to the situation in 
which the field is parallel to the c-axis. In the case of having magnetic field applied parallel to 
the graphene planes nearly no MR is observed. This fact speaks for a huge anisotropy of the 
metallic or superconducting regions and it suggests that those are within the interfaces found 
by TEM, see Fig. |2l We note that the field-induced metal-insulator transition (MIT) found in 
bulk HOPG samples, see Section 1421 vanishes for samples thinner than 50 nm, which also 
indicates that certain regions parallel to graphene planes are related to the origin of this MIT 
tBarzola-Ouiquia et aU.2008,) . 



Will-b8-s8t-by-IN-TECH 



■31 00 



■ 




i 



CD 



1 



10 



100 1000 
Sample Thickness t(nm) 



10000 



Fig. 5. Resistivity at 4 K vs. sample thickness of different samples obtained from the same 
HOPG bulk sample. 

2.4 Thickness dependence of the resistivity 

Figure |5] shows the resistivity at 4 K of different samples vs. their thickness. It is clearly 
seen that the resistivity decreases increasing the sample thickness. The average change in 
resistivity between ~ 10 nm to 17 fim thick samples is about two orders of magnitude, far 
beyond geometrical errors. A possible explanation for the observed behavior is that the larger 
the thickness the larger the amount of interfaces in the sample, see Fig.|2] As we described 
in Section 12.31 there is a clear change in the temperature and magnetic field dependence 
decreasing the sample thickness. It appears unlikely, however, that randomly distributed 
point-like lattice defects can be the reason for the observed behavior Also the interpretation 
provided by Zhang et al. ( 2005) that the decrease of mobility (i.e. an increase in the resistivity 
at constant carrier density) decreasing sample thickness provides an evidence for boundary 
scattering is surely not the correct explanation for the observed behavior, taking into account 
the fact that one graphene layer shows finite and large mobility. Note that thin graphite flakes 
show ballistic transport with huge mobilities, see Section |3] The overall results suggest the 
existence of a kind of thickness-threshold around ~ 50 nm for multigraphene samples of few 
micrometers size, obtained from HOPG bulk graphite indicating that the interfaces contribute 
substantially and in parallel to the graphene layers. 

In the past the influence of the internal interfaces in graphite bulk samples was either 
completely neglected or the scientific community was not aware of their existence. On the 
other hand it is well known that grain boimdaries in various semiconductors can lead to the 
formation of quasi-two dimensional carrier systems confined in the boundaries. As early 
examples in literature we refer to the quasi-t wo dimensional electron gas system that was 
found at the inversion layers of n-G e bicrystals jUchida et allll983l:IVul & Zavaritzkaval.ll979 ) 
or in p-InSb ([Herrmann et alj.ll984h as well as in Hgi _ TCdxTe grain boundaries jLudwig et al , 
Il992h . It is important to note that in Ge-bicrystals lUchida et al.l Il98^) found actually the 
quantum Hall effect (QHE) at T ^ 4.2 K and at magnetic fields above 6 T. The density 
of carriers at the interface was estimated to be ~ 5 x 10^^ cm^^. Therefore, we note that 
the usually reported carrier concentrations for graphite are not intrinsic of ideal graphite. 
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as sh own by lArndt et alJ l |2009h and also by iDusari. Barzola-Ouiquia. Esquinazi & Garcial 
furthermore, we speculate that the QHE behavior observed in bulk graphite samples 
Kempa et alJ.l2006l:lKopelevich et aD.l2003h comes from internal interfaces. 

3. Searching for the intrinsic transport properties of the graphene layers inside 
graphite 

3.1 Background 

In this Section we will briefly discuss the background of the ballistic transport and its 
experimental observation in graphene layers inside the graphite structure. The wave nature 
of the electrons plays an important role when the sample dimensions are comparable with the 
wavelength of the electrons. This turns to be possible in graphite samples because the density 
of carriers is very small, increasing therefore the Fermi wavelength Af . On the other hand, as 
we will see below, a carrier within the graphene layers inside graphite can transit micrometers 
through the sample without collisions. Therefore the carriers in the graphene layers within an 
ideal graphite structure have the unusual property of having large Af as well as large mean 
free path £. 

The transport phenomenon in the ballistic regime is best described by the Landauer approach. 
Consider a narrow constriction connected through two wide contact leads across which a 
voltage is applied. Let the electrochemical potentials of two regions be /(j and /(2- The net 
current flowing through the device can be taken as I = {2e^/h) I]T;(/^i — 112) where T,- is 
the transmissio n probability of carriers and the factor 2 comes f rom the spin degeneracy in 
each sub-band ( lButtikeii[l986l:lLandauerl,ll957l:lTsu & Esakilll973h . The effective conductance 
through a narrow constriction is given by G = 1 / — }i2) = {le^ /h)Y,Ti- When the 
constriction dimensions are comparable to Af and much shorter than the mean free path £, one 
can see steps in the conductance when the constriction width is varied. These steps are nearly 
integer multiples of {2e^) /h. Conductance quantization can be also achieved by varying the 
potential energy. As the constriction width is reduced or the gate voltage is changed in a 
determined direction, the number of propagating modes at a given energy decreases, i.e. the 
sub-bands are cut off one by one, and therefore the conductance decreases. A new mode 
appears in the conductance when the constriction width increases by Af / 2. 
Nowadays, it is feasible to fabricate devices that show ballistic transport. However, it 
is difficult and time consuming to fabricate these devices because of the required small 
dimensions when common metals or semiconductors are used. In the past there were some 
experiments reveling conductance quantization at low temperatures. The first experimental 
observation s of conductance qua ntiza tion in two dimensio nal electron gases (2DEG) were 
reported by Ivan Wees et all jl988^ and I Wharam et aP h988|) . They studied the point contact 
conductance in GaAs-AlGaAs heterostructures as a fimction of a negative gate voltage at 
low temperatures. The conductance showed clear plateaus at integer multiples of le^ /h as 
the width increases by an amount of Af /2. In both cases the authors could show that the 
transport is completely ballistic and the conductance is determined by the number of occupied 
sub-bands, independently of the channel length. 

3.2 Ballistic transport in graphite 

For decades, the transport properties of bulk graphite were interpreted using 
Boltzmann-D rude approach assuming diffusive ohmic behavior for the conduction electrons 
teeliyl Il98lh . Within this approach one has four temperature-dependent free parameters, 
namely two mobilities and two carrier densities. The determination of these parameters 
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using this model implies fitting transport data with at least four free parameters. To obtain 
accurate values for these parameters it is necessary, however, to go beyond this model. 
When the mean free path and Fermi wavelength are of the order of sample size, one is not 
allowed to use the Boltzmaim-Drude transport theory to determine the electrical resistance. 
As discussed before if the size of the system is of the order or smaller than the carrier mean 
free path, ballistic regime enters in which carriers can move through the system without 
experiencing any scattering. Usually in metals this takes place in the nanometer range. 
However, graphite is extraordinary because its mean free path is o f the order of microns 
jPusari. Barzola-Ouiquia. Esquinazi & GarciaLl201ll;lGarcia etal 1 1200^ . 
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Fig. 6. Scanning electron microscope pictures of the two multigraphene samples, A (left) and 
B (right) with the Pd electrodes. Sample A has a constriction of 4f/m and sample B of 1.5 fim. 
The scale bar in sample A indicates 10 pim and in sample B 5 pim. The dashed red lines 
denotes the border of the multigraphene samples. 

There is a transparent method to obtain all the main transport parameter s from exper iments 
based on constrictions prepared in the graphite or multigraphene samples jCarcia et al.l, l200S). 
When the mean free path is lar ger than the co nstriction width ( W), the transport of the carriers 
shifts from ohmic to ballistic l |Sharvinl.ll965h . By measuring the longitudinal resistance as a 
function of the con striction width, one can determine the mean free path, carrier density and 
Fermi wavelength jPusari. Barzola-Ouiquia. Esquinaz i & Garci4 l201ll;[Garcia et"ar lllooi). 
Consider the geometry depicted in Fig. 1 of iGarcia et al. (20qB), or the samples in Fig. (6] 
for a constriction in a quasi-two dimensional sample. Already Maxwell pointed out that a 
constricting circular orifice of diameter W in a thin, nontransparent screen of size L, produces 
a spreading resistance that is equal to the p{T)/W if the ratio £/W << 1 (p is the resistivity 
of the material). However, when this ratio increases there are two corrections to the Maxwell 
spreading resistance: (i ) the Ohmic value is corrected by a factor of the order of unity as 
pointed bv lWexletlil966h , and ( i i) more irnportant, a dominant b allistic term appears. This was 
observed by lGarcia et alj j2007| ); Knudsen ' 1934 ); SharvinI jl965l) and the value of the resistance 
R in three dimensions reads Icarcia et alil2007l;IWexlerlll966h : 



^3D 



ApiT)£ , j(K)piT) 



3A 



W 



(1) 



where A = 7tW^/4 is the area of the hole or constriction and 7(k) is a smooth varying 
geometrical function, i.e. 7(k = W/£) ~ 1 — 0.33/ cosh(O.lK) = 0.67 . . . 1 for A: = O...00 
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jWexleil Il966h . In Eq. (l) the first and the second terms of the right-hand side (rhs) 
correspond with the Knudsen-Sharvin and Ohmic resistances, respectively. The spreading 
Ohmic resistance in three dimensions can be estimated within a factor 2/7T off from the exact 
Maxwell's solution assuming a hemisphere ir\ which the electric field E(r) = /3D(r)p(T). The 
radius r is taken at the constriction middle point and J^oi^) is the current density equal to 
the total current I divided by half of a sphere, i.e. hoi^) = U (27rr^), assuming that due to 
symmetry the current is radial far away from the constriction. 

From a similar calculation but in the two dimensional cas e, appropriate for graphite due to 
the weak coupling between the graphene planes, we have dGarcla et al 

R (T\-A^\PIT) I , 2p(r)7(K)/n(0/W) p{T)L 

R2d{T) - a—e{T) + a + (2) 

The first term at the right-hand side (rhs) of Eq. (|2) corresponds to the ballistic 
Knudsen-Sharvin resistance; the second, logarithmic term to the ohmic, spreading resistance 
in two-dimensions; here n is the total sample width and t its thickness. The logarithmic 
dependence on the constriction width of this ohmic, diffusive contribution is due to the 
quasi-two dimensionality of the transport in graphite and supports the assumption of weakly 
interacting graphene layers inside the sample. The last term is due to the ohmic resistance 
of the constriction tube itself. From previous works it was clear that the position and 
shape of the voltage el ectrodes affects the outcome of the experiment in mesoscopic devices 
jMcLennan et al.l, Il99ll) . Therefore the constant a was introduced, which takes care of the 
influence of the sample shape, the topology, and the location of the electrodes in the sample. 
For the usual electrode positions through the whole sample width as shown in Fig. \6\ 
a = 1. The validity of Eq. (|2), especially the logarithmic dependen ce of the ohmic par t, for 
HOPG as well as for multigraphene samples has been verified by iGarcla etal and 
iDusari. Barzola-Ouiquia. Esquinazi & Garcial (|201lh . 

In the following, we review some of our experimental results for two multigraphene 
samples, A and B, with different geometry and resistivi ty. The sample details a s 
well as their preparation and fabrication were de scribed by iBarzola-Ouiquia et alj ( |200^ ; 
iDusari, Barzola-Ouiquia, Esquinazi & Garcial boill) . The constrictions in the middle of the 
samples, see Fig.|6l were prepared with the focused-ion beam of a dual-beam microscope. It is 
important to note that we avoided the modification of the crystalline structure of the samples 
due to the ion beam spread by protecting them wit h a ^ 300 nm thick negative-e lectron beam 
resist (AR-N 7500), a method successfully tested by lBarzoIa-Ouiquia et alJ l|2010f) . 
To obtain the mean free path without further adjustable parameter we measured the resistance 
R as a function of the constriction width W and use Eq. (|2). Figure [7| shows the results for 
samples A and B at two temperatures. The results show that for W < 1 iim the ballistic 
contribution (dashed lines in Figs, ^a) and (b)) overwhelms the ohmic ones indicating that 
the mean-free path should be of this order. Having only £ as free parameter Eq. (|2) can be used 
to fit the behavior of R vs. W for sample A. From the theoretical fits one obtains £ = 1.2 jim 
and 0.8 iim at 60 K and 250 K, respe ctively. The Fermi wavelength per graphene layer can be 
calculated using dGarcla et al. lEo^: 

1 2np{T)e{T)e^ 

^' = v, h ' 



where = 0.335 nm is the distance between graphene planes in the graphite Bernal stacking 
configuration. For sample A we obtain then Af = 0.5(0.8) ± 0.25 fim at 250 K (60 K). 



12 



Will-ba-s8t-by-IN-TECH 



Despite of the go od agreement obtained for sample A and for HOPG bulk samples 
iCarcia et all |2008|) , Eq. ^ suffers from an important limitation since it carmot describe 
correctly ballistic transport phenomenon in which the wave nature of the electrons plays a 
crucial role, i.e. in samples where Af > W. As discussed before, in this case the ballistic 
contribution to the resistance is better described by the quantization of the transverse electron 
momentum in the constricted region. In this case the value of the resistance is given by 
the inverse of a sum of an energy-dependent and trans verse wave vectors ij„ -dependent 
transmission probabilities T„, where n = 0, ±1, ±2, . . .Nc jStauber et all 1200^ . These values 
correspond to the one dimensional electric sub-bands and Nc is the largest integer smaller than 
2W / Af . In this case the increase in resistance is expected to show an oscillatory behavior as a 
function of W or Af jGarcia & Escapal 1989 : Snvman et alll2003) as observed experimentally 
in bismut h nanowires jCosta-Kramer et all 19971 ) as well as in GaAs devices ivan Wees et all 
[1988; Wh aram et al"lll988h . 

Note that the obtained mean free path for sample A is smaller than the distance between the 
electrodes. The larger the sample, larger is the probability to have defective regio ns with larger 
carri er concentration and smaller mean free path within the voltage electrodes jArndt et all 
l2009l) . In order to increase the probability to observe the expected quantization phenomenon 
in multigraphene samples, it is necessary to have a mean free path larger than the sample size 
in order to be completely in the ballistic regime. Therefore, we repeated the experiment with 
sample B that shows lower resistivity and with a smaller voltage-electrode distance, see Fig.[6] 
Figures [Ztc) and (d) show the measured resistance normalized by its value at a constriction 
W = 3 fim for sample B. The normalization is necessary because in this way we pay attention 
to the huge relative increase of R decreasing W and we need neither the absolute value of 
p nor of a to compare the data with theory. We realize that for sample B Eq. (2j does not 
describe accurately the experimental data even assuming the largest possible mean-free path 
equal to the voltage-electrode distance of ~ 2.7 ^m. The data can be reasonably well fitted 
dividing the ballistic term in Eq. (2j by the function trunc(2W/AF)Af /2W, which generates 
steps whenever the constriction width W ~ i\f / 2 with i an integer. From the fits we obtain 
the parameters Af = 1.0(1.5) ^m and £ = 2.2(2.7) ± 0.3 }im at 300(10) K. Using other values 
of £, for example E = 1.3 ;(m, see Fig.[7td), the function does not fit the data indicating indeed 
that the carriers behave ballistically between the voltage electrodes, leaving actually Ap the 
only free parameter. 

The ballistic analytical function we used resembles the theoretical results with similar steps 
obtained by lSnyman et alj j2008l) where the conductance vs. W was calculated numerically 
for a single layer graphene with an electrostatically potential landscape that resembles a 
constriction. An analytical average value or envelope of this stepped function is obtained 
replacing the truncation function by exp(— Af /2W), see Figs.[7l!c) and (d). This exponential 
function represents the impossibility of an electron to propagate in the constriction when 
W < Af/2. The important result obtained for sample B is the huge increase of the 
res istance for W < 2 ^m ind icating clearly a larger E than the one obtained in sample A, 
see lDusari, Barzola-Ouiquia, Esquinaz i & Garcial feoill) for further details. 
The temperature dependence of R(T, W) can be used now to obtain Af(T) and the mobility 
per graphene layer, this last given by fi(T) = {e/h)Af(T)£(T). Since the density of carriers per 
graphene layer can be calculated from n = In/ we show in Fig. [5] the mobility vs. carrier 
density for the two samples and for a bulk HOPG sample and compare them with data from 
literature for suspended single layer graphene. From these results we clearly recognize the 
much larger mobility and smaller density of carrier for the graphene layers inside graphite, 
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Fig. 7. Measured resistance as a function of the constriction width W at (a) 60 K and (b) 250 K 
for sample A. The point with the largest W corresponds to the virgin sample without a 
constriction. The different lines correspond to the ballistic contribution (dashed line) and 
ohmic second (dotted line) and third (dashed-dotted line) terms in Eq. l|2). The continuous 
line represents the addition of the three contributions. For (a) the continuous line is 
calculated with £ = 1.2 /im and for (b) (. = 0.8 f(m. (c) and (d): Normalized resistance for 
sample B vs. constriction width W at 10 and 300 K. Note the semi logarithmic scale. The line 
with steps is obtained dividing the ballistic term in Eq. l|2) by (AF/2W)trunc(2W/AF) with 
the parameters £ = 2.7 fim and Af = 1.5 jim. The dashed line follows Eq. (2). The dotted line 
is obtained multiplying the ballistic term in Eq. {2) by the exponential function exp(AF/2W). 
(d) The same as in (c) but the continuous hne was obtained with the parameters Ay = 1.0 fim 
and (. = 2.2 ^m. The dashed stepped function is obtained using the same Ap but with a 
smaller £ = 1.3 }im. Adapted from Dusari, Barzola-Ouiauia, Esauinazi & Garcia l2 01lh . 

supporting the view that the graphene layers within graphite are of better quality and with a 
smaller carrier density than single layer graphene. Obviously neither sample A nor B nor the 
HOPG sample are free from defects and therefore we expect that the obtained values might 
still be improved in ideal, defect-free graphite structures. 
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Fig. 8. Mobility vs. carrier density obtained for sam ples A and B (after 
Dusari. Barzola-Quiquia. Esquinazi & Garcial l|201lh ) and for HOPG (dashed line, after 
Garcia et al.l bOOsI) '). Note that the data points are not obtained at a constant temperature but 
the temperature is changed to get a change in the carrier density and mobility. The points are 
obtained between 300 K and 10 K for sample B and 300 K and 60 K for sample A. For the 
HOPG sample (dashed line) the data run from 5 K at the largest mobility values to 300 K. Th e 
straig ht line corresponds to the data of a suspended graphene sample at 20 K from lDu et al.l 
<200a) . 

4. The influence of interfaces inside bulk graphite samples: Hints for the 
existence of high-temperature superconductivity 

4.1 General remarks 

Ideal graphite is a layered material where each of the layers (graphene) consists on a 
honeycomb lattice of carbon atoms. Graphene layers are stacked together with Van der 
Waals forces much smaller than the covalent ones between the carbon atoms within a single 
layer. This weak coupling between graphene layers within the graphite structure can be 
"translated" in a huge anisotropy in the resistivity pel pa ^ 10^ at low temperatures. In high 
quality samples this leads to a quasi-two dimensionality of the tr ansport because most of 
it occ urs within the graphene layers inside the graphite structure jKopelevich & Esquinazil 
|2007[) . However, real graphite cannot be considered as an ideal unif orm stack of ideal 
graphene layers, as we disc ussed in Section [2l The role o f defects jArndt et all. l2009l; 
iBarzola-Ouiauia & E sQuinazi l2010l;lBarzola-Ouiquia et ail|2003) as well as internal interfaces 
Barzola-Ouiquia et aL, ,2008i) have a crucial effect in the electrical properties. Clear evidence 
on the non uniformity of real graphite structures is given by several microscopy techniques. 
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see Section 12] In this Section we are interested in particular on the influence of the internal 
interfaces on the transport. 

4.2 Bulk graphite 

Starting with the bulk material, iKempa et alj l|2000l) studied the magnetoresistance in the 
in-plane direction in HOPG bulk samples where a magnetic-field-driven transition from 
metallic- to semiconducting-type behavior of the bas al-plane resistance pg was foimd. Thi s 
was reproduced in later publications by other groups | |Du et alll2005l:lTokumoto et alJ.l2004h . 
Note, however, that depending on the thickness of the graphite samples and their quality, 
metallic or semiconducting behavio r in graphite can be observed if no magnetic field is 
applied ^Barzola-Ouiquia et ID. l2008h . There is a clear suppression of the metallic-like phase 
(in case it appears) b y a magnetic field applied perpendicular to the graphene planes (i.e. 
parallel to the c-axis) <Kempa et Later work found a field-induced metal-insulator 

transition (MIT) also in the c-axis resistivity of graphite pc jKempa et al.l. |2002|) . The 
metallic-like behavior of pc has been related to the one found in the longitudinal resistivity 
due to a conduction-path mixing mechanism, i.e. the conduction path of the carriers along the 
c-axis is in part short circuited by lattice defects and/ or impurities. 

If non-percolative superconducting grains are located at interfaces running mostly parallel to 
the graphene planes, the largest increase in the resistance with magnetic field applied parallel 
to the c-axis is expected in the temperature region where either the resistance decreases 
decreasing temperature or it levels off at zero field applied. This is expected since at a 
temperature below the maximum or leveling off one speculates that the coupling between 
superconducting grains starts to be observable. Kopelevich et al. (1999) showed a sensitive 
change in the temperature dependence of the resistivity in bulk HOP G samples imder an 
applied magnetic field parallel to the c-axis. Additional studies done by lKempa et al.l i2003|) 
show the absence of a MIT in graphite if the magnetic field is applied parallel to the graphene 
planes, a fact that indicates that its origin is related to regions running parallel to the graphene 
layers. 

4.3 Anomalous hysteresis loops and quantum resonances in the magnetoresistance 

In percolative homogenous superconducting samples, one can ascribe superconductivity by 
observing a screening of the external magnetic field (Meissner effect) below a critical field 
and /or by measuring the drop of resistance to practically zero below a critical temperature. 
Inhomogeneous granular superconducting materials (in which the superconducting phase 
could be found in small parts of the whole sample) require other way of testing, as resistance 
does not drop necessa rily to zero and the Meissner effect is probably too small to be measured. 
lEsquinazi et al.l feOOSi) found anomalous hysteresis loops in the MR similar to those observed 
in granular supercond uct ors with Josephson-coupled grains, see Fig.[9t a). 
lEsquinazi et al.l ^20081) and lDusari. Barzola-Ouiquia & Esquinazil i201l|) show that the minima 
in the MR are located in the same quadrant field from which one start s the field sw eeping 
as it was reported in conventional high Tc granular superconductors ili et all Il993h . This 
can be explained based on a two-level critical-state model where pinned fluxons exist inside 
the Josephson-coupled superconducting grains but also between them. These last ones are 
usually much less pinned and therefore can strongly influence the MR behavior, especially at 
low enough fields. The fact that the minimum of the MR is rather flat can be interpreted as 
due to non-uniform superconducting grain size. 
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In addition, it must be stressed that when the field is parallel to the graphene planes no 
effect is observed, suggesting that the superconducting regions or patches are rather parallel 
to the graphene layers. Quantum resonances in the MR were also observed, see Fig. IHb), 
and ascribed to Andreev's reflections between localized superconducting patches connected 
by semiconducting r egions in which Cooper pairs may flow relatively large distances 
jEsquinazi et From the amplitude and period of the oscillations in fi eld, the distance 

between superconducting granular domains was estimated to be < 1 }im jEsquinazi et in. 
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Fig. 9. (a) Weak hysteresis in the resistance of a thin graphite sample vs. magnetic field, 
similar to that observed in superconductor s with josephson-coup led grains near zero field, 
measured in a multigraphene sample (after lEsquinazi et alj feOOSl) '). (b) Oscillations in the 
resistance of a multigraphene sample as function of the magnetic field at 2 K after subtraction 
of a linear field backg round. The continuous line corresponds to a Fourier fit. Adapted from 
lEsquinazi etaTI <200d) . 



4.4 Behavior of the magnetoresistance of multigraphene samples with micro-constrictions 

The novel method of using thin mesoscopic samples with micro-constrictions not only 
allows us to obtain basic parameters like the mean free path and carrier density, see 
Section |3l but provides also the possibility to increase the sensitivity of the voltage 
measurement to the superconducting regions that may be at the constriction region. We 
mean that the expected superconducting patches inside the constriction area might be better 
"detected" by usual voltage measurernent a t two contacts on either side of the constriction. 
jPusari, Barzola-Ouiquia & Esquinazil boill ) verified this effect and foimd an increase in 
the range of temperature where the previously reported fesquinazi et 

iUIIool) 

anomalous 

hysteresis loops can be observed, see Fig.lTOl 

All the experimental results indicate that, in case that the superconducting regions exist, 
they should be localized mainly at the interfaces between the crystalline regions of slightly 
different orientation observed by TEM, see Fig. |2l At these interfaces the density of 
carriers should be high er iough to achieve high criti cal temperatures, provided the quasi-two 
dimensionality remains jGarcia & Esquinazil, l2009l) . The sample can be modeled as the 
sum of superconducting and normal in series and in parallel circuit paths. A coherent 
superconducting state at low enough temperatures within the superconducting patches is 
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Fig. 10. Relative difference between the two field branches of the hysteresis loop in the MR 
(see Fig.|9la)) as a function of temperature for two different muitigraphene samples at a 
constant field o f 1.66 mT. Black squares (■) correspond to the data presented by 
Esquinazi et al. ilOOSi); (*, O, *) correspond to the results presented by 
Dusari. Barzola-Ouiquia & Esquinazil l|201lh for a muitigraphene sample without and with 
two constrictions. The increase in the hysteresis using the constriction method is clearly seen. 
The samples without constrictions show granular superconducting-like hysteresis loops up 
to 10 K only, while samples with the constrictions show this behavior up to 50 K in the case of 
4 }tm constriction width. 

expected as well as in between some of them. Thermal fluctuations influence the coherent 
superconducting state in some parts of the sample and no zero resistance is achieved. To 
summarize the results reported so far in tens of nanometers samples as well as the first 
indications in HOPG bulk samples, two fundamental aspects must be remarked: first, 
metallic-like behavior of the resistivity in HOPG bulk is not intrinsic and not related to the 
scattering of conduction electrons with phonons, and second, if superconductivity plays a 
role in the anomalous properties observed in HOPG samples then the q uasi-two dimensional 
interfaces may contain the superconducting regions dGarcia et 



4.5 Transport measurements of TEM lamellae 

As explained before, the indications for superconductivity in graphite get clearer when the 
sample dimensions are reduced to some extent. On the other hand, if the sample dimensions 
are further reduced, the intrinsic properties of the graphene layers inside graphite can be 
investigated. In order to further investigate the role of interfaces a different way of preparing 
samples from bulk graphite seems necessary. The main idea is to prepare small and narrow 
graphite samples (within the graphene planes) in order to confine the path of the input current 
to go through some of the possible superconducting interfaces running mostly parallel to the 
graphene planes. Further details of the process as well as the lamella dimensions can be seen 
in Fig.[lT] Several lamellae have been studied and the results are presented in the following 
lines. 
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Fig. 11. Sketch of the lamellae preparation process. Focused Ion Beam (FIB) is used to cut 
bulk HOPG (up left side: reddish areas represent the area to be cut with Gallium ion beam 
and this last is represented by the red arrows). Previous to this cut, the surface of HOPG is 
protected with WO thick layer deposited using the electron beam. Lamellae of size 
~ 7 ^m xl5 jim x200 nm (this last dimension in the graphene plane direction) are prepared 
as shown in the upper right draw. The low-TEM picture (bottom right) shows the interfaces 
between the crystalline regions (one of them marked in red). Blue lines drawn in the upper 
pictures provide an idea where the quasi-two dimensional superconducting patches are 
supposed to be located. The bottom left picture represents how the situation at one of the 
interfaces might be. The blue areas simulate the superconducting patches that might be 
weakly, Josephson coupled by low-carrier density graphene regions. 

In Fig.[l2]we show the change of the voltage vs. temperature measured at two positions in two 
different lamellae, samples 1 and 2, both prepared using the same procedure. The estimated 
critical temperatures for these samples are: 25 K and 175 K, respectively, although we should 
clarify that these temperatures not necessarily are the intrinsic critical temperatures of the 
patches but probably those related to the Josephson coupling between them. The differences 
observed in the values of Tc reflect the inhomogeneous structure expected for the graphite 
interfaces. Sample 1 results can be considered as the "ideal ones" because within resolution 
zero values in the voltage are reached for low enough temperatures and small input currents, 
see Fig. ll2r a). In this sample, the used electrical contacts for the voltage measurements were 
directly at or near some superconducting paths that short circuited below certain temperature 
and input current. Sample 2 shows similar results, qualitatively speaking, but the range of 
critical values in temperature and current are different. In this sample the contacts are not 
placed directly on the superconducting regions and therefore normal regions of the samples 
contribute to the whole measured voltage. It can be seen in Fig. \V2\ h) that in this case the 
transition is sharper and a change of two orders in magnitude in the measured voltage is 
achieved. Due to the van der Pauw contact configuration used, negative values (respect to 
the input current direction) of the voltage below the transition temperatures are observed for 
several lamellae, instead of zero voltage. This can be explained using a Wheatstone bridge 
circuit as it will be shortly described later in this Section. We found a strong dependence of 
the measured voltage on the input current. Above a certain current neither negative nor the 
drop in voltage are observed anymore. 
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(a) Sample 1 (b) Sample 2 

Fig. 12. Voltage vs. temperature measured for two different lamellae. Note the clear change 
in the temperature dependence with input current. A drop in the measured voltage is 
observed for all lamellae when the input current is small enough. 

Current- Voltage characteristic curves (CVC) are the next important piece of evidence to be 
mentioned. In case of having Josephson coupling the expected CVC should show a nearly 
zero voltage value in the measured voltage below a certain value of current and at low 
enough temperatures. At currents much above a critical one an ohmic linear behavior should 
be recovered. Another important feature that particularly at low temperature is found in 
Josephson-coupled systems is the irreversibility in the CVC. As shown in Fig. I13f a) this is 
achieved in the case of the "ideal" sample 1. The CVC results for other two lamellae become a 
bit more complex but compatible with the previous results shown in Fig. [T2l b). Figures lTSf b) 
and (c) correspond to the CVC for samples 2 and 3, respectively. Negative values in the 
differential resistance are measured below certain values of current and temperature and 
above them a linear ohmic behavior is recovered. The simplest way to explain these results 
is by using a non-homogeneous current-voltage circuit composed by four different resistors 
Rj,i = 1 . . .4 within a simple Wheatstone bridge model as for example: V{I) oc I[R^R^ — 
R2R3) / {{Ri + R2){R3 + R4)). Each of the Rj represents an effective resistance within a path 
inside the lamella. In case one of the resistors includes Josephson-coupled superconducting 
regions with zero voltage value at low enough temperatur es and currents, we can simulat e 
the measured behavior using f or this R, the CVC given by I Ambegaokar & Halperinl ( Il969h . 
lAmbegaokar & HalperinI jl969t) described the voltage vs. current behavior within the DC 
Josephson effect including thermal noise. 

As it can be seen in Figs. [T3l b) and (c) the theoretical calculations (solids lines) - using one 
free parameter (the critical Tose phson current) to calculate the corresponding following 
lAmbegaokar & Halperinl il969l) and fixing two or three arbitrary current-independent resistor 
values for the rest of the Wheatstone circuit - appear to fit well the experimental data. 
Obviously, pure Josephson-coupled resistors are rather difficult to be observed experimentally 
because the measured voltage is the sum of contributions coming from the superconducting 
as well as the semiconducting regions where the first ones are embedded. The case of sample 1 
is observed in ~ 20% of the measured lamellae. 

As shown in Section 14.3! the role of the magnetic field in the granular superconductivity 
behavior is definitively of huge relevance. When superconducting patches are present 
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Fig. 13. Current-voltage characteristic curves at different fixed temperatures for three 
different lamellae. The continu ous lines that fit the data were obtained using the 
lAmbegaokar & HalperinI il969l) CVC for two R,'s in (b) and one R,- in (c) within a 
Wheatstone bridge model, assuming different critical Josephson currents at different 
temperatures. No magnetic field was applied and the earth field was not shielded. 

between semiconducting normal regions (we assume that they run mostly parallel to the 
graphene layers) a remarkable effect must appear when a magnetic field is applied in the c-axis 
direction and less influence should be observed if it is applied parallel to the graphene planes. 
The results in Fig.[14]correspond to the situation where the field is applied parallel to the c-axis 
of the lamellae. The in-plane case is not shown here because no effect has been observed, as 
expected. Figure iMf a) shows the temperature behavior of the measured voltage with and 
without field. The shift to lower temperatures of the transition after applying 0.75 T is clearly 
shown for one of the samples. In the case of other sample a 0.1 T field radically changes 
the previously observed drop in voltage and a semiconducting-like behavior appears. This 
behavior resembles the MIT found already in HOPG and discussed in Section l4!2l Figure lTU b) 
shows the CVC at 4 K and 25 K with and without field. One can recognize that the effect of 
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(a) (b) 

Fig. 14. (a) Voltage dependence on temperature for two different lamellae. Upper picture: 
The transition denoted by the sharp decrease in the voltage below a certain temperature 
(black squares) is shifted to lower values when 0.75 T is applied perpendicular to the planes 
(red dots). Bottom picture: Similar to the other picture but for a different lamella. The zero 
field (black squares) transition is completely suppressed after applying 0.1 T (red dots), (b) 
Current-voltage characteristic curves for a different lamella, with (red dots) and without field 
(black squares) at 4 K (upper graph) and 25 K (bottom graph). 

the field is to reduce the superconducting effect, i.e. the CVC tends to a linear, positive voltage 
behavior. 

iHeersche et alj j2007l) showed that Cooper pairs can flow in a single graphene layer relatively 
long distances between two superconducting contacts. The authors showed also the existence 
of Josephson currents. These results support the main idea proposed in this Section 
concerning the weakly Jose phson coupled s uperconducting regioris throu gh a graphene layer. 
In addition to the work of iHeersche et alj i2007h . iHagymasi et alj (|2010h developed a theory 
for the Josephson critical current in ballistic Superconductor-Graphene-Superconductor 
structures. We use their model to calculate the reduced temperature dependent critical 
Josephson current. In Fig. [15] the reader can see the solid line obtained from the model 
proposed and the experimental data taken from the experimental work described in this 
Section. The model fits well the data. In the near future further investigations concerning 
the effect of the magnetic field in these lamellae are needed. 
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Fig. 15. Normalized Josephson critical current Ic / kiO) as a fimction of the normalized 
temperature T/Tc obtained for different samples. The values obtained for lc{0) range from 
0.9 to 7 jiA and for Tc ~ 23 k to 175 K for different s amples. The solid line is the theoretical 
curve taken from the work of iHagymasi et all i2010l ) without free parameters. 

5. Conclusion 

The main messages of the work described in this chapter are: (1) Several of the transport 
properties measured in large graphite samples are not intrinsic of graphite but are influenced 
in a large extent by interfaces and other defects and impurities. Taking into accotmt that one 
defect or impurity in a graphene layer inside graphite can contribute with one carrier, it is clear 
that the carrier densities are not intrinsic but are due to defects and to a large extent probably 
to highly conducting interfaces, internal as well as at the sample boundaries. We note also 
that, as the large experimental evidence in semiconductors indicates, even the quantum Hall 
effect observed in graphite is actually not intrinsic of ideal graphite but can be attributed 
to the highly conducting quasi-two dimensional electron systems at the interfaces. Clearly, 
the details of this effect will depend on the sample dimensions and quality, providing us a 
way to understand experimental discrepancies found in literature. Our work indicates that 
the intrinsic mobility of the carriers in the graphene layers inside graphite are extraordinarily 
large /; > 10^ cm^/Vs and these carriers show ballistic transport in the micrometer range even 
at room temperature. 

(2) The first observation of superconductivity in doped graphite goes back to 1965 when i t 
was observed in the potassium graphite intercalated compound C8K jHannay et al.l. [l96^ . 
A considerable amoimt of stu dies has reported this phenomenon in intercalated graphit e 
compounds or doped graphite temery et al"ll2005l:lKopelevich et aLl,l2004l:IWeller et al.l,|200l) , 
however the superconducting properties of pure graphite are still imder discussion. In this 
chapter we have given a possible clarification of this topic. By reducing the dimensions 
of the samples, the intrinsic properties of graphite have been deeper investigated. The 
overall idea coming out from all the experimental results can be resumed as follows: bulk 
graphite samples cannot be considered as a uniform electronic system, but it must be 
considered as a semiconducti ng matrix with metallic as well as superconducting domains 
in it. We note that already ILu et al.l l|2006h realized the coexistence of insulating and 
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conducting behaviors in bulk graphite surfaces. From this newly acquired knowledge the 
observed behaviors appear compatible with the existence of non-percolative partially highly 
conducting, partially superconducting regions coupled by Josephson-coupling through 
graphene planes. The interfaces observed by TEM may have enough carrier density to trigger 
quasi-two dimensional superconductivity. 

Theoretical work that deals with superconductivity in graphite as well as in graphene has 
been published in recent years. For example, p-type superc onductivity has been predicted to 
occur in inhomogeneous regions of the graphite structure iGonzalez et al.l, 200 1[) or d— wave 
high -Tc superconductivity based on resonance valence bonds jBlack -Schaffer & Donia"chl. 
I2OO7I) . Following a BCS approach in two dimensions critical temperatures ~ 60 K 
have been obtained if the density of conduction electrons per graphene plane increases to 
n 10^^ cm~^, a densi ty that might be induced by defects and /or hydrogen ad-atoms 
jGarcia & Esquinazill2009l ). Further predictions for sup erconductivity in graphene support the 
premise that n > 10^^ cm^^ in order to reach Tc > 1 K feopnin & SoninLl2008l;IUchoa & Netol 
12007 '). The interfaces observed by TEM might then be the regions where enough carrier 
density exists to trigger quasi-two dimensional supercond uctivity. In contrast to the basically 
3D superconductivity in intercalated graphitic compounds jCsanyi et 3111200^ we expe ct that 
supe rconductivity at quasi-2D graphite interfaces as well as at doped surfaces I Ha n et aP. 
120101 ) may exist at much higher temperatures , partially because of the r ole of the high-energy 
phonons in the 2D graphite structure itself jCarcia & EsquinazH |2009|) . Room temperature 
superconductivity with n d + id pairing symmetry ha s been predicted to occur in doped 
graphene with a carrier concentration n > 10^^ cm^^ jPathak et all I2OI0I) . We note that if 
we take into account the density of interfaces then a measured carrier density in bulk graphite 
samples of 5 x 10^^ cm^^ would mean an effective density ~ 2 x 10^^ cm~^ at the int erfaces. 
And last but not least we refer to a recent theoretical work from lKopnin et all j201ll) where 
the authors emphasize that a topological protected flat band in semimetals may promote 
superconductivity at very high temperatures. 
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